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Filling factor v = 1 bilayer electron systems in the quantum Hall regime have an excitonic- 
condensate superfluid ground state when the layer separation d is less than a critical value d c . On 
a quantum Hall plateau current injected and removed through one of the two layers drives a dissi- 
pationless edge current that carries parallel currents, and a dissipationless bulk supercurrent that 
carries opposing currents in the two layers. In this paper we discuss the theory of finite supercurrent 
bilayer states, both in the presence and in the absence of symmetry breaking inter-layer hybridiza- 
tion. Solutions to the microscopic mean-field equations exist at all condensate phase winding rates 
for zero and sufficiently weak hybridization strengths. We find, however, that collective instabilities 
occur when the supercurrent exceeds a critical value determined primarily by a competition between 
direct and exchange inter-layer Coulomb interactions. The critical current is estimated using a local 
stability criterion and varies as (d c — d) 1//2 when d approaches d c from below. For large inter-layer 
hybridization, we find that the critical current is limited by a soliton instability of microscopic origin. 

PACS number: 73.40. Hm, 73. 20.Dx 



I. INTRODUCTION 

In bilayer quantum Hall systems, broken symmetry 
ground states [1,2] that have spontaneous interlayer 
phase coherence were predicted some time ago [3-8] . This 
broken symmetry state is expected to be most robust 
near total Landau level filling factor v = 1 and occurs 
only when interactions between electron in opposite lay- 
ers are comparable in strength to interactions between 
electrons in the same layer. The putative existence of this 
broken symmetry was used several years ago to explain 
[9] the observation of a surprisingly strong dependence 
[10] of the bilayer system v = 1 charged excitation gap 
on in-plane magnetic field strength. 

Spontaneous coherence between electrons in different 
energy bands is an old topic in condensed matter physics, 
although it has not yet been convincingly demonstrated 
outside of the quantum Hall regime. For example, it 
has long [11] been realized that spontaneous coherence 
at zero magnetic field is a possibility when overlapping 
or nearby bands have opposite quasiparticle energy vs. 
wavevector curvatures, the conduction and valence bands 
in a semiconductor or a semimetal in particular. Recent 
studies of optically generated electron-hole plasmas in 
semiconductors, do indeed hint [12] at the expected col- 
lective behavior. In a separate materials system, the dis- 
covery of weak ferromagnetism in lightly-doped divalent 
hexaborides [13], which are ferromagnetic despite the ab- 
sence of partially filled d- or /- orbitals, led Zhitomirsky 
et al. [14] to propose recently that spontaneous coher- 
ence between conduction and valence bands could be the 
mechanism responsible for their ferromagnetism. Spon- 
taneous coherence between different hyperfine states in 
Bose-Einstein condensates of magnetically confined 87 Rb 
atoms, manifested by Rabi oscillations between the two- 



components, is also closely related to the quantum Hall 
bilayer phenomena studied here. [15] 

Spontaneous coherence states are most commonly de- 
scribed using the language of semiconductor physics in 
which a particle-hole transformation is performed for the 
valence band; spontaneous coherence between the bands 
then maps to electron-hole pair condensation, something 
that is closely analogous to Cooper pair condensation in 
a superconductor. These ordered states can be equiva- 
lently described as a pseudo-spin- 1/2 quantum ferromag- 
nets with the two states (top or bottom layer label in a 
quantum Hall bilayer system, band-index in a bulk semi- 
conductor, or alkali atom species label in an ultra-cold 
atom system) defining an ordered spinor [9] . The ordered 
state has superfluid properties for "staggered" currents 
that flow in opposite directions in the two bands, exci- 
tonic superfluidity in the language of semiconductors. It 
was suggested more than twenty years ago [16] that spon- 
taneous coherence could occur between two-dimensional 
(2D) conduction and valence bands localized in separate 
quantum wells and that [17,18], at least in mean-field 
theory, the conditions required for condensation were 
more likely to be met when the 2D systems experienced 
a strong perpendicular magnetic field. It has not always 
been recognized, however, that because of the dispersion- 
less Landau bands that occur in a strong magnetic field, 
there is no difference [4] in this regime between spon- 
taneous coherence between a conduction band and a va- 
lence bands and spontaneous coherence between two con- 
duction bands. The spontaneous coherence that occurs 
in bilayer quantum Hall systems is, in fact, precisely that 
originally anticipated by early theoretical work. [17,18]. 
Advances in our understanding of quantum Hall physics 
have, however, given us a deeper appreciation of the lim- 
itations of the mean-field theory approach used in the 
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older work and of the physics behind its partial success. 
Superflow of the electron-hole pair condensate is perhaps 
the most characteristic property of the bilayer quantum 
Hall broken symmetry states, which we refer to here as 
quantum Hall superfluids. In this paper we discuss pre- 
dictions for the maximum sustainable (counter-flowing) 
staggered supercurrents supported by these states that 
follow from microscopic mean-field theory for the ideal 
case in which the two-dimensional electron layers are 
completely free from disorder. 




FIG. 1. A schematic illustration of a bilayer quantum Hall 
bar. For v = 1 and sufficiently small interlayer spacing d this 
system exhibits spontaneous interlayer phase coherence and 
excitonic superfluidity that supports counter-flowing (stag- 
gered) supercurrent J s . The spiral is intended to indicate a 
uniform interlayer phase gradient that leads to constant J a 
along x. 

The work presented here complements an earlier field- 
theoretic investigation by Kyriakidis and one us [19], who 
studied the decay of staggered supercurrents in bilayers 
quantum Hall superfluids by thermally activated phase 
slips and calculated the staggered I-V characteristics and 
critical current as function of interlayer gate voltage. The 
main technical difference between the two studies is that 
in the work of Ref. [19] it was assumed that the phase 
slip instability of a super-current carrying state is at long 
wavelengths, an assumption that turns out to be valid 
only for screened short-range interactions. As we explic- 
itly show here, for realistic (long-range) Coulomb inter- 
actions, the instability is at a finite wavevector. Fur- 
thermore, we argue that even with short-range electron- 
electron interactions, when the interlayer spacing d is 
near the critical value d c (determined by the transition 
out of interlayer coherent quantum Hall state), the in- 
stability is also at short microscopic scales of order the 
magnetic length. This finding limits the range of validity 
of the field theoretic approach of Ref. [19] to the case of 
very closely spaced layers, which cannot yet be realized 
experimentally. 

Recently experimental advances by Eiscnstcin and col- 
laborators [20-23] have revealed several dramatic effects 
that are only partially understood but are believed to be 
due to collective transport effects in quantum Hall super- 
fluids. It appears, however, that the samples studied in 
current experiments have only short-range order because 
of disorder. Their progress, which we believe will open 
up a set of subtle new questions about non-equilibrium 
properties of superfluid-likc ordered states, provides pow- 



erful motivation for the present work. The key capabil- 
ity which allows collective particle-hole transport effects 
to be probed by electrical experiments is the possibility 
of making separate electrical contact to two-dimensional 
electron gas layers with a separation d w 20nm, small 
enough that is to be in the range where spontaneous 
coherence occurs. The strong zero-bias peaks they see 
in interlayer tunneling conductance [21] studies partially 
confirm predictions [24,8,25] of Josephson-like effects, al- 
though these experiments are not fully understood and 
two of us has argued elsewhere that the analogy with the 
dc Josephson effect is incomplete. [26,27]. The evolution 
of the tunneling peak when the in-plane field [22] is var- 
ied shows evidence of the predicted [3,7,25,28] linearly 
dispersing Goldstone collective mode associated with ex- 
citonic superfluidity. Most directly related to the present 
theoretical paper, is a very recent experiment [23], which 
studies transport properties of a quantum Hall superfluid 
for the case in which current (integrated across the sam- 
ple) flows through only one of the two layers. Closely 
related transport phenomena have also been discussed 
theoretically [29] for the case of thin film ferromagnets. 
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FIG. 2. Schematic non-equilibrium quasiparticle popula- 
tions for the quantum Hall effect in a Hall bar geometry for 
an ideal disorder-free quantum Hall bilayer superfluid. The 
quasiparticle states are labeled by a guiding center coordinate 
which varies over the range from one edge of the sample to 
the other and is proportional to canonical momentum along 
the Hall bar. The quasiparticle current is equally divided 
between the two layers and a net current is maintained by 
having different chemical potentials on the two edges. The 
quasiparticles have spontaneous phase coherence and carry 
edge currents that are evenly partitioned between the two 
layers. In order to satisfy the drag experiment constraint 
that no net current to flow in one of the layers, the conden- 
sate must carry an electron-hole supercurrent equal in magni- 
tude to the quasiparticle current. Unlike the charge current, 
the electron-hole supercurrent in the ideal case will flow uni- 
formly through the bulk of the system. The mean-field state 
of a quantum Hall superfluid pairs electrons and holes that 
have different momenta in the current carrying direction, or 
equivalently, different cyclotron-motion guiding centers. This 
property is indicated schematically by the slanted lines that 
connect different guiding centers. 
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The quantum Hall effect and superfluidity share the 
unusual property of transport without dissipation. In su- 
pcrfluids, dissipationless transport is possible because the 
quasiparticles are in equilibrium with a current-carrying 
condensate. In the quantum Hall effect, as illustrated 
schematically in Fig. 2, the quasiparticles arc not in equi- 
librium. Instead the gap for charged excitations in the 
bulk implies that low-energy quasiparticles are localized 
at the sample edges, allowing a net current to be car- 
ried through the system [30] without dissipation by main- 
taining a Hall voltage difference between isolated subsys- 
tems on opposite edges. In the mean-field theory [31] 
of a quantum Hall supcrfluid, the occupied quasiparticle 
state wavefunctions are coherent linear combinations of 
orbitals localized in separate layers. For equal density 
in the bilayers, the current they carry is divided equally 
between the layers and any voltage probe that couples to 
the quasiparticle system will measure the same value in 
either layer, leading to large drag voltages [32]. Since the 
quasiparticles carry equal current in the two layers, the 
only way in which it is possible to have no net current in 
one of the layers, is to have it canceled by spontaneously 
generated staggered supercurrent, as illustrated in Fig. 3. 
In fact, as discussed in Ref. [33,34], since both staggered 
supercurrent and the uniform Hall current are dissipa- 
tionless, at v — 1 an ideal quantum Hall supcrfluid state 
should exhibit, respectively, vanishing and quantized lon- 
gitudinal and Hall drag-resistivities. 




FIG. 3. Schematic for an integer quantum Hall superfluid 
bilayer in the drag geometry of Kellogg, et al., Ref.23, where 
electrical current Ji is only allowed to flow in the top layer. 
To satisfy the boundary condition of vanishing current in the 
bottom layer, J2 = 0, the system spontaneously develops a 
dissipationless staggered current J s . 

Current experiments are not fully in accord with this 
simple picture; the voltage drops in the current carrying 
direction are not zero and not identical in the two layers, 
and the Hall voltages are not perfectly quantized as one 
would expect in a v = 1 quantum Hall bilayer superfluid 
at T = 0. Most likely the discrepancy is due to disor- 
der which leads [35] to mobile quasiparticle excitations 
in the bulk of the two-dimensional electron system and 
to vortex flow and consequent phase slips in the super- 
fluid order parameter field. Long length scale inhomo- 
geneity in the samples that leads to phase separation be- 
tween ordered and disordered regions can also be invoked 



to explain [36] many aspects of current experimental re- 
sults. The calculations we report in this paper do not 
account for disorder and do not attempt to explain quan- 
titative aspects of current experiments. Instead we focus 
on the properties of ideal systems in which no disorder is 
present, with the expectation that the ideal situation will 
be approached more closely in the future. We expect the 
collective instability studied here to control longitudinal 
dissipation and deviation from Hall-drag quantization in 
clean quantum Hall bilayer supcrfluids. 

In Section II of this paper we briefly summarize the 
mean-field-theory of current-carrying states in quantum 
Hall supcrfluids, allowing for the possibility of explicit 
symmetry breaking inter-layer hybridization terms in the 
microscopic Hamiltonian and for weak links in the two- 
dimensional supcrfluid that are created by gate voltages. 
We define the critical current as the maximum current at 
which stable local minima in the energy-functional of the 
quantum Hall superfluid exist; formal expressions for the 
stability matrix at mean-field energy-functional extrema 
are also given in Section II. Some of these formal mean- 
field-theory results have been discussed in another con- 
text in earlier work, [37] but are included here for com- 
pleteness. In Section III we apply these formal results 
to the case of a uniform quantum Hall superfluid with 
no inter-layer hybridization. In this case we are able to 
derive physically transparent expressions for the critical 
current. An important prediction that follows from this 
analysis is the way the critical current vanishes as the 
phase boundary that separates the quantum Hall super- 
fluid from normal ground states of the 2D electron sys- 
tems is approached. Section IV deals with the potentially 
interesting case [38] in which weak links are intention- 
ally induced in these two-dimensional electron systems 
by applying voltages to external gates. We find that even 
relatively weak disturbances lead to drastically reduced 
critical currents. This partially explains the sensitivity 
of experimental samples to disorder, that will undoubt- 
edly introduce such weak links in the superfluid bilayer. 
In this case we are able to derive a Josephson-like rela- 
tionship between the order parameter phase jump across 
the weak link and the staggered supercurrent that flows 
across it. 

An important difference between excitonic supcrflu- 
ids and both superconductors and fermion-pair super- 
fluids is the inevitable presence in the Hamiltonian of 
explicit staggered gauge-symmetry breaking terms that 
allow electrons to tunnel between layers and hybridize 
quasiparticle states even when there are no interactions. 
These terms can be extremely weak but are never strictly 
zero. In Section V we show that staggered current causes 
a soliton lattice to form at a characteristic current den- 
sity. For very strong tunneling, we find that the maxi- 
mum supercurrent drops to zero. We conclude in Section 
VI with a summary of our results and a discussion of the 
role of quenched disorder and thermally activated vor- 
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tices and phase slips that is informed by our numerical 
results for the microscopic properties of ideal systems. 
We conclude that non-linear transport in the drag exper- 
iment geometry used by Eisenstein and collaborators [23] , 
studied as a function of carrier density, might be able to 
identify a crossover in dominant dissipation mechanism 
that is closely associated with the ideal system critical 
currents that are estimated in this paper. 



II. THE SUPERFLUID STATE ENERGY 
FUNCTIONAL: MEAN-FIELD EQUATIONS AND 
STABILITY ANALYSIS 

All of the calculations summarized in this article 
are based on a microscopic Hartree-Fock approximation 
(HFA) for the energy functional of quantum Hall stag- 
gered superfluid states and are analogous to ones made 
using the BCS theory for superconductors. We will con- 
sider only situations in which the order parameter is a 
function of a single (X) spatial coordinate, making it 
convenient to choose a Landau gauge in which the Lan- 
dau level orbitals are localized in this direction. Be- 
cause we ignore edge effects, it will prove more con- 
venient to choose orbitals localized along the direction 
of current (X) flow rather than the more conventional 



choice perpendicular to the direction of current flow il- 
lustrated in Fig. 2. The energy functional is Emf[9, <p] = 
(*|T + V ext + V ee \$>) where T is the interlayer tunnel- 
ing term, V ex t is the external potential term, and V ee is 
the Coulomb interaction in the microscopic Hamiltonian. 
The variational wavefunction used to define this energy 
functional is 

|*) = [] (cos[9(X)/2]c XT + sin[0(X)/2]e*W4 B ) |0), 
x 

(1) 



where X is a Landau gauge guiding center label, T 
(top) and B (bottom) are layer indices, and 9 and ip de- 
fine the coherence factors of the BCS-like single-Slater- 
determinant many-particle wavefunction. If the phase 
coherent state is regarded as a pseudospin fcrromagnet, 
the angles 9{X) and <p{X) are polar and azimuthal an- 
gles that specify the direction of the pseudospin magne- 
tization at position X . This is a special case of a class of 
variational wavefunctions in which virtual charge density 
fluctuations are not permitted. The energy functional 
can be written explicitly [40,41,37] in terms of Hartree 
and exchange microscopic two-particle interaction matrix 
elements: 
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where As as is the splitting between symmetric and antisymmetric orbitals induced by hybridization, V g (X) is the 
external gate voltage averaged over guiding center orbital X, and H, Fs, and Fjj are the Hartree matrix element, the 
same- layer exchange matrix element, and the different-layer exchange matrix element. If we neglect the finite widths 
of the quantum wells the interaction matrix elements are given by: 



dq 27re 2 (l - e - qd ) 
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(3b) 



where Vc(q x , %) = 2ire 2 /(eq) and 27re 2 exp(— qd)/(eq) for 
C = S and C — D respectively. Note that the exchange 
integral drops rapidly with orbit center separation, while 
the electrostatic integral falls only as X~ 2 at large X, 
corresponding to the interaction between dipole moment 
lines created by charge transfer between the layers. In 
these equations £ = yJhc/eB is the magnetic length. The 
Fourier transform of these functions will figure promi- 
nently in analyzing the superfluid properties of uniform 
bilayers. The energy functional Emf[0(X),(p(X)] can 
be elevated to a quantum theory by adding to it a Berry 
phase contribution to the Lagrangian [9,39]. 



We choose the convention 
dX 



H(p) 



-J 



2kP 



cnp(i P X)H(X) 



(4) 



(and correspondingly for the exchange integrals) so that 
these quantities have units of energy When finite- 
thickness effects are neglected, we obtain the following 
explicit expressions for the Coulomb interaction case: 
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where Jo(x) is the Bessel function. The same layer ex- 
change function Fs(p) differs only by the absence of the 
er' rd factor. F D (p) is plotted in Fig. 4. 

A. Mean-field Equations 

Mean-field states of quantum Hall supcrfluids are ex- 
trema of the energy functional (2) and satisfy mean-field 
equations obtained by setting 8Emf/5^p{X) = 0, and 
SEmf/S0(X) = 0. These Euler-Lagrange equations for 
the variables 9 and <p are coupled, and highly non-linear. 



taxnp(X) = 



Our primary interest here is in current carrying states, 
and we choose a to be the distance over which the conden- 
sate phase increases by 2ir; our mean-field states are ones 
in which phase slips form a lattice. Note that sin(ip(X)), 
cos((p(X)) and r](X) are periodic. It follows that the to- 
tal phase change of the order parameter on going from 
one end of the system to the other is 2nL/a = 2nN w 

These equations must be solved self-consistently to lo- 
cate mean-field states. They can be understood most 
simply in the language of pseudo-fcrromagnetism, in 
which they simply state that at each position X, the ori- 
entation of the quasiparticle pseudo-spinor is along the 
direction of the total pseudo-spin effective field, including 
direct and exchange mean-field interaction and external 
potential contributions felt by the quasiparticles at that 
point. The quasiparticles of the superfluid are in equilib- 
rium with the condensate at extrema of the mean-field 
energy functional. 

B. General Expression for the Condensate Current 

The expression we use for the condensate staggered su- 
percurrent follows from arguments presented in greater 
detail in previous work on quantum Hall supcrfluids [9] 
and implicitly in work on excitonic superfluidity. The 
current operator in each layer can be expressed in terms 
of the derivative of the Hamiltonian with respect to 
the vector potential in that layer. The vector poten- 
tials in the two layers can be varied by independent 
gauge transformations for the two layers and physical 
properties can depend on the phase difference between 
the two layers only through the gauge-invariant quantity 
hVip + (e/c)(AT — Ab) where At and Ab are the vector 
potentials in top and bottom layers and ip(X) is the dif- 



The equations are, in addition, strongly non-local be- 
cause of the long range of the Coulomb interaction term. 
For V g (X) and As as equal to zero, the global minimum 
of the energy functional is 0(X) — ir/2 and <p(X) equal to 
an arbitrary constant. Since we will be interested primar- 
ily in solutions with 9(X) near 7r/2, it is useful to define 
rj(X) = it/2 — 0(X). For our investigation we will be in- 
terested only in extrema that are periodic functions of X. 
Denoting the period by a and using the thermodynamic 
limit property J2x ~~ * I dX/(2n£ 2 ), the mean-field 
equations can be written as: 



(6a) 
(6b) 



ference in wavefunction phase between top and bottom 
layers that appears in the HF variational wavefunction, 
Eq.l. Since the operator for the spatially averaged cur- 
rent in a layer is proportional to the derivative of the 
Hamiltonian with respect to the vector potential in that 
layer and we have chosen the same gauge for the vector 
potential for each layer, it follows that the spatially av- 
eraged condensate current is proportional to the rate of 
change of condensate energy with N w : 

j _ 1 dEuF _ a 2 deMF 
s nh L y dN w nh da 

where £mf is the energy per cross-sectional area. This 
is the expression that we use below to extract conden- 
sate staggered currents from solutions of the mean-field 
equations. In Eq. 7 J s is a number staggered current den- 
sity, i.e. it corresponds to an electrical current density 
equal to — eJ s (X) in the top layer and eJ s (X) in the bot- 
tom layer. When tunneling between the layers is allowed, 
there is a circulating current within each period of the 
phase slip lattice in addition to this spatially averaged 
current. The circulating current is conveniently calcu- 
lated from charge conservation which relates the current 
flowing between the layers to the divergence of the con- 
densate supercurrent. The experimental meaning of the 
condensate current evaluated in this way is discussed at 
greater length in Section VI. 

C. Stability Analysis 

We restrict ourselves to the case of zero bias voltage 
for which the mean field solutions will have rj(X) = 0. 
A mean-field solution is a local minimum of the energy 
functional provided that the stability matrices 
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are both positive definite. The explicit expressions for 
the stability matrices evaluated at the mean field solu- 
tions rj{X) = 0, and ip = <p(X) arc as follows: 



2H(X - X') - F S (X - X')} 



Asas cos(^(X)) + -L - X ") C0S M x ) - ¥>(*")) 
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K<p<p{X, X') 



-T^FniX - X') cos^(X) - <p{X')) 



>X.X' 



A SAS cos(ppO) + ±- ]T F D (X - X") cos(^(X) - ip{X")) 



(10) 



As we show in the next section, for Agyis = 0, the 
mean-field solutions do not break translational invari- 
ance, and the stability matrices can be diagonalized ex- 
plicitly by taking advantage of the translational symme- 
try property. The spectrum of the stability matrices can 
then be related to the Fourier transformed direct and 
exchange interaction matrix elements as we show below. 
In the more general case of finite Asas, it is necessary 
to evaluate their spectrum numerically, which we do in 
Sec.V by discretizing the set of allowed guiding centers. 



III. CRITICAL CURRENT OF UNIFORM 
QUANTUM HALL SUPERFLUIDS (A S as = 0) 

For Ag^s = and V g (X) equal to a constant, there are 
solutions to the mean-field equations for which r\ (bilayer 
charge imbalance) is spatially uniform and ip(X) varies at 
a constant rate. The largest possible supercurrents flow 
for balanced bilayers, i.e. for V g (X) = 0. The family of 
solutions with rj = and ip = QX corresponds to states 
in which electron-hole pairs have condensed into a state 
with total momentum Qx. Unlike the case of the zero- 
field fermion pair condensates, solutions to the mean field 
equations can be found for any pairing momentum Q, no 
matter how large. Inserting these solutions in Eq.(2), we 
find that the condensate-dependent part of the energy 
per unit area for pairing at wavevector Q is Fd(Q)/8tt£ 2 . 
The condensate current then follows from Eq.(7) 
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-1 2dF D (Q) 



2 d,£MF 

h dQ ~ 4nhe 2 dQ 
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The dependence of J s on pairing wavevector Q is illus- 
trated in Fig. 4 for a range of bias voltages. The natural 
unit of charge current density in the quantum Hall regime 
is J = 2e 3 /he£ 2 w 82^A/^mB[Tesla]. 
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FIG. 4. J s vs. pairing wavevector Q for Asas = V g 







and a series of layer separation d values. J a is in units of 
J /e ee (Jo = 41 x 10 _6 Amp/^mB[Tesla].) Our sta- 

bility analysis demonstrates that the mean-field state is al- 
ways unstable for pairing wavevectors larger than Q* , the 
pairing wavevector at which J s is maximized. For quantum 
Hall superfluids, however, the critical current is almost al- 
ways limited by unstable amplitude fluctuations of the elec- 
tron-hole pair condensate that are indicated by negative val- 
ues of Kzz(q). For example for d — £, the maximum super- 
current mean-field state (J c = 0.0058 Jo) occurs at Q£ = 1.6, 
but amplitude fluctuations are unstable for Q£ > 1.2 as we 
discuss below. The inset shows Fd(Q) vs. Q for d = £. 
Fd(Q) ~ constant — Q 2 at small Q and ~ l/Q at large Q so 
that J s is proportional to Q and Q~ 2 at large and small Q 
respectively. 



The main effect of a uniform bias voltage is to reduce 
the scale for supercurrent values. We should also remark 
that these supercurrent values and all numerical results 
discussed in this paper will be influenced somewhat by 
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corrections that account for the finite thickness of the 
two-dimensional layers. These corrections could be in- 
corporated into our calculations without any difficulty, 
but we choose to ignore them here mainly for the sake 
of simplicity. As we discuss later, other more fundamen- 
tal and difficult issues arise when we attempt to compare 
this analysis with experiments in real systems. 

The Hamiltonian Hf for Gaussian fluctuations <j>(X) = 
<p — QX and m z (X) around these uniform superfluid 
(<po = QX, m z = 0) mean-field states can be expressed in 
terms of the stability matrices discussed previously. In 
the uniform current case, K lpip (X,X') and K ZZ (X, X') 
depend only on X — X' and the matrices can be diago- 
nalized by Fourier transformations. 

We find that 



m z (X). (p(q) = Y.x<t > { x ) ex P{- i P x )l N 4> wh ere iV = 
L x L y /(2irl 2 ) is the total number of Landau gauge or- 
bitals in the system. This choice is made so that <p(q) 
is dimensionless and (p(q = 0) is the spatial average of 
4>(X).) The first two terms in Eq.(13b) represent the cost 
in electrostatic energy and the gain in exchange energy 
that accompanies fluctuations in the bilayer charge bal- 
ance. Both stability kernels K vv and K zz must be posi- 
tive definite for the superfluid state to be stable against 
small fluctuations. From Eq.(13b), we sec that for q — > 0, 



K ir . 



Afiq 



2 dJ B 
dQ 



which is positive only on the increasing 



portion of the J S (Q) curve; at larger pairing wavevectors 
phase separation into regimes with larger and smaller 
average phase winding rates is always energetically pre- 
ferred over the homogeneous state. 
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FIG. 5. K zz (q) and K vv (q) vs. q for d = £, and Asas = 0. 
K vv {q = 0) = even at finite pairing momentum Q / 0, be- 
cause the energy is invariant under a spatially constant change 
in the phase tp. K zz (q) is nonanalytic at q = 0, decreasing 
in proportion to \q\d, because of long-range Coulomb interac- 
tions which favor non-uniform exciton density variations. The 
decrease of K zz (q) for qd small is rapid for large layer sepa- 
rations d, and eventually drives the critical current to zero, 
destabilizing the superfluid state. 

(In these equations we have used the following 
conventions for the Fourier transforms of (f>(X) and 



FIG. 6. Critical supercurrent density J c (d) (main panel) 
and (inset) pairing wavevector Q c (d) of the critical current vs. 
layer separation d, defining the QH Superfluid phase bound- 
ary in pairing wavevector-layer separation space. The dashed 
lines in the main figure and the inset showing critical current, 
J„ and maximum pairing wavevectors, Q, obtained by ap- 
plying the K vv (condensate phase) order parameter stability 
criterion is shown by these results to be less stringent; the 
solid line is determined by the K zz stability criterion. For 
d w d c , J c oc (d c — d) 1 / 2 , while for d — > J c oc d. Cur- 
rent experimental samples have layer separations close to the 
critical value d c . Because the K zz (q) instability is at a fi- 
nite wavevector g», we expect that for small d the resulting 
current carrying state is a modulated QH superfluid, i.e., a 
QH supersolid. For larger values of d, the instability is at a 
large finite wavevector and we expect a first order transition 
between uniform superfluid and disordered ground states. 

For q — > and Q = 0, K vv — > p s q 2 where p s is phase 
stiffness (superfluid density) parameter of the superfluid. 
K vv (q) appears to always be positive when it is positive 
for q — > 0; the K vv stability criterion Q < Q* (with 
Q* the pairing wavevector at which J S (Q) is maximum) 
does not limit the critical current of quantum Hall su- 
perfluid to values below the maximum value at Q* that 
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occurs in mean-field solutions. In contrast, K zz {q) al- 
ways has its minimum at a finite value of q. This dif- 
fers from the field-theoretic treatment of Ref. [19], where 
K zz {q) was taken to be a monotonically increasing func- 
tion of q. This assumption is not valid for unscreened 
(long-range) Coulomb interactions. The last term in the 
expression for K zz (q), which represents the loss in con- 
densation energy when the condensate electron-hole den- 
sity is varied from its optimal value, is essential for the 
stability of quantum Hall supcrfluids at typical values of 
d/£. As illustrated in the inset of Fig. 4, this stabiliz- 
ing term decreases in magnitude at finite Q because of 
the reduced condensation energy associated with finite 
Q pairing. It is this reduction in condensation energy 
that limits the critical staggered supercurrent in our the- 
ory. It is also worth remarking that K zz {q) is negative 
for any finite value of Q for d = 0. It is only at d = 
that the ground state is given exactly by the Halperin 
[1] (1, 1, 1) two-component quantum Hall fluid wavefunc- 
tion. It follows that the (1, 1, 1) state is not a superfiuid, 
a point that we have made [42] previously. Because of 
this property, field-theoretical analyses of quantum Hall 
superfiuid properties that start from (1,1,1) variational 
wavefunctions should be regarded with some caution, in 
our view. 

In Fig. 6 we have plotted the layer separation depen- 
dence of the ideal bilayer critical current J c that follows 
from this stability analysis. The critical current vanishes 
as d — > because of the vanishing difference between in- 
tralayer and interlayer interactions that is necessary to 
provide a barrier to phase slip nucleation. In the lan- 
guage of pseudo-ferromagnetism [9,50], the correspond- 
ing observation is that easy-plane anisotropy is required 
for spiral state metastability. The staggered critical su- 
percurrent also vanishes at large layer separation because 
the stiffness against electron-hole (bilayer charge imbal- 
ance) density waves vanishes as we discuss below. This 
figure shows that the critical current is always limited 
by the K zz {q) stability criterion. The maximum critical 
current occurs for d ~ QAt and is ~ 0.008 of our current- 
density unit Jo. It follows from this analysis that under 
typical experimental circumstances the critical current 
of an ideal system should be ~ 0.33/zA/zm _1 . This stag- 
gered critical current value is several orders of magnitude 
larger than the value used in recent drag experiments 
[23]. Although existing drag experiments provide evi- 
dence of collective transport, they explicitly demonstrate 
dissipation in the pseudospin current channel because the 
measured longitudinal electric fields in the two layers are 
different. This dissipation is likely due to the flow of vor- 
tices induced by thermal fluctuations and disorder that 
we do not account for here. It would nevertheless be in- 
teresting to experimentally explore the possibility that 
a large increase in the difference in longitudinal elec- 
tric fields would occur relatively sharply at larger current 
densities, reflecting a change in the dominant dissipation 



mechanism when the ideal critical current addressed in 
this paper is exceeded. 

The K vv (q), and K zz (q) stiffnesses in the fluctua- 
tion Hamiltonian are inversely related to static response 
functions. Since the phase of the condensate and the 
electron-hole pair density are canonically conjugate co- 
ordinates these two stiffnesses are also simply related to 
the system's quantized collective excitation energies. The 
semiclassical picture of these modes emerges from the 
Landau-Lifshits equations of motion of the pseudospin 
ferromagnet, which can be derived by adding a Berry 
phase term to the above fluctuation Hamiltonian to ob- 
tain the system's Lagrangian [9] . 
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FIG. 7. The dispersion relation of the quantum Hall super- 
fluid collective modes is plotted for d = £, and V g = Asas = 
for values of the pairing wavevector Q which vary between 
Q — (the top curve) and Q = 0.7-T 1 . At Q = Q c , where 
the roton minima vanishes, the superfiuid mean-field state is 
unstable. 

The superfiuid collective modes have energy e{q) = 
2ir£ 2 K zz (q)K ipip (q) , generalizing earlier results 
[18,3,28,9] to the case of a finite supercurrent state. 
For all pairing wavevectors, the collective modes have a 
linear dispersion of a superfiuid at long wavelengths. At 
shorter length scales (large q) e(q) also has a minimum, 
that is reminiscent of the roton minimum in superfiuid 
Helium. The source of this characteristic minimum in 
the collective mode dispersion e(q) is the Coulomb inter- 
action, which leads to the minimum in K zz {q) displayed 
in Fig. 5. The collective mode dispersion is altered when 
a supercurrent is flowing primarily because of the re- 
duction in the condensation energy and K zz (q) with 
increasing Q. When the minimum is already close to 
zero for Q = 0, the critical supercurrent is small. Since 
Fd{Q) is a quadratic function of Q for small Q and the 
minimum in K zz (q,Q = 0) approaches zero linearly as 
d — > d c , it follows that Q c and the critical staggered 
supercurrent vanish like (d c — d) 1 / 2 for d — > d c . 
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The physics that limits the ideal critical current of 
quantum Hall supcrfluids is similar to that encapsulated 
by Landau's expression for the critical current of super- 
fluid Helium, in which he argued that the superfiuid ve- 
locity cannot exceed J™ ax = v™ ax n = min (e*(q) /hq)n, 
where e*(q) is the elementary excitation dispersion re- 
lation when no supercurrcnt is flowing, i.e. at Q = 0, 
and the minimum above is over all wavevectors q. In 
both superfiuid Helium and quantum Hall superfiuid 
cases, Landau's argument approximates the more fun- 
damental requirement that the order parameter of the 
current-carrying state be at a local minimum of the 
model's energy functional. The remark that we make 
here applies mutatis mutandis to conventional super- 
fluids. We first note that for small Q, K zz (q,Q) rts 
K zz (q,0) — p s Q 2 so that the maximum value of Q is 
[min(K zz (q, 0))/p s ] 1//2 . One additional approximation is 
required to obtain the same result from the Landau cri- 
terion; we must assume that K vv (q,Q = 0) ~ p s q 2 for 
q up to the value at which s*(q)/q is minimized, actu- 
ally a reasonable approximation for quantum Hall su- 
pcrfluids. Using n = (2n£ 2 )~ 1 , the density of a full 
Landau level, we see that the superfiuid velocity is re- 
lated to Q by v s = p s Q{2-k£ 2 )/Ti. The Landau crite- 
rion v 2 < min[e*(<7) 2 /<7 2 ] w (2ir£ 2 ) 2 mm[K zz (q, 0)p s ] then 
leads to the same limit on Q as the K zz (q, Q) > stabil- 
ity criterion, that we use here to obtain the QH superfiuid 
phase boundary illustrated in Figs. 6, 7. 



IV. WEAK LINKS IN QUANTUM HALL 
SUPERFLUIDS (A S as = 0) 

One element of the physics of quantum Hall supcrflu- 
ids that might potentially be interesting for future ex- 
perimental study, is the possibility [38] of creating and 
tuning weak links using gate voltages. Weak links are 
locations in the superfiuid at which the local superfiuid 
stiffness is reduced and phase changes are therefore con- 
centrated in a narrow region (link). In quantum Hall 
bilayer superfluids, one such weak-link geometry is illus- 
trated in Fig. 8, where the local stiffness is reduced in situ 
by a gate voltage that locally moves the bilayer system 
away from balanced condition of equal density in the two 
layers. 

In this section we discuss the supercurrent properties 
of quantum Hall superfluids with such top gate-voltage 
induced weak links. We represent the electrostatic po- 
tential from a narrow gate separated from the epitaxially 
grown bilayer system by a cap layer (Fig. 8), somewhat 
crudely, by an external potential that has a narrow Gaus- 
sian variation in the x direction, but is invariant in the 
y direction, as shown in the inset of Fig. 9. The width 
of the Gaussian should be given approximately by the 
minimum of the gate width and the thickness of the cap 
layer. [47] 
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FIG. 8. Schematic illustration of a quantum Hall superfiuid 
weak-link induced by a top gate. The pseudospin stiffness is 
reduced locally by a gate voltage that displaces the bilayer 
system away from the balanced condition of equal density in 
the two layers. The increased phase winding rate at the link, 
required by current conservation, is illustrated. 



In order to evaluate the Josephson-like relationship be- 
tween the phase change across a weak link and the su- 
perfiuid current, we have solved the mean-field equations 
for a periodic system system with one weak-link and a 2ir 
phase change per period a. In this way we can vary the 
superfiuid current by changing the period a. A typical 
solution is shown in Fig. 9. 
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FIG. 9. <fi(X) and m z (X) are shown for a bilayer with 
Asas = and d — £. These results are for a weak link po- 
tential with a width equal to I and a barrier height 1.5e 2 /e£. 
Aip is measured with respect to its value at one edge of the 
periodic cell and the weak link is located at the center of the 
well. 



Since for vanishing interlayer tunneling, As as = 0, 
considered here the dc staggered supercurrent is spatially 
uniform, current continuity demands dip/dX to be in- 
versely related to the local superfiuid density when the 
phase gradient is small. In particular, dip/dX must be 
constant to an excellent approximation away from the 
weak link. We will refer to this constant as Q below. 
After solving the mean-field equations for a series of Q 
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values, the current density can be evaluated either by ex- 
ploiting its relationship to the a dependence of the total 
energy or by combining the phase winding rate Q away 
from the weak link with the uniform system J S {Q) ex- 
pression. The phase change across the weak link, Aip, is 
extracted by extrapolating the constant curve from one 
side of the weak link to the other and finding the addi- 
tional shift that is present because of the weak link, as 
illustrated Fig. 9. 

A typical result for J s vs. Aip obtained in this way is 
summarized in Fig. 10. 
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ical current decreases rapidly with both the width of the 
barrier potential and its strength. Because sample in- 
homogeneities will undoubtedly introduce a distribution 
of weak links into the superfluid bilayer, one conclusion 
that follows immediately from these calculations is that 
the critical current will likely sharply reduced by disorder 
in typical current samples. 
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FIG. 11. Critical staggered supercurrent density depen- 
dence on the maximum of the bias potential y™ ax responsible 
for the weak link, plotted for two width values. 



FIG. 10. Supercurrent density vs. weak-link phase slip an- 
gle Aip for weak links of width I and heights V g max = 0.75 
and V™ ax = 1.50e 2 /e^. The maximum supercurrents in this 
case are approximately five times less than in the absence of 
a weak link. 

For small currents (large a), J s increases linearly with 
Aip, as in a conventional Josephson junction. Notice that 
J s reaches its maximum for Aip close (but not exactly 
equal) to tt/2 as in the Ambegaokar-Baratoff [44] theory 
of superconductors linked by a tunneling Hamiltonian. 
The microscopic physics of the coupling is quite different 
in the quantum Hall superfluid, however. For a weak link 
that is translationally invariant along y [47] (the geom- 
etry illustrated in Fig. 8) the guiding center X remains 
a good quantum number in the absence of disorder and 
consequently the quasiparticles do not tunnel across [47] . 
Instead, the coupling between the uniform superfluids on 
opposite sides of such weak link is due to Coulomb ex- 
change interactions, both within and near to the weak 
link. As in the case of a uniform superfluid, configura- 
tions of the order parameter field that occur beyond the 
maximum of this function are not stable, and the stag- 
gered critical current is given by the maxima of these 
curves. 

In Fig. 11 we plot critical currents as a function of the 
strength of the bias potential that produces the weak link 
for two different values of the width of the link. The crit- 



We close this section by noting that to obtain the weak- 
link critical current J Cl we have relied solely on the van- 
ishing K w criterion, and have ignored the K zz stability 
matrix. Although we have not checked this explicitly, 
because a weak link substantially reduces the phase stiff- 
ness, leading to a five-fold reduction in </*, we believe it is 
this K w criterion that determines the weak-link critical 
current, giving J c = J*. 

V. ROLE OF INTERBAND HYBRIDIZATION IN 
QUANTUM HALL SUPERFLUIDS: THE A S as + 
CASE 

In the presence of interlayer tunneling As as, the su- 
perfluid energy functional has soliton lattice extrema 
which break translational symmetry but, as we comment 
later, still carry a finite supercurrent. It is not in general 
possible to address the metastability of these solutions 
analytically, although a qualitative understanding can be 
achieved using scaling arguments [19]. 

A. K vv instabilities of charge-balanced bilayers 

To find members of this family of energy functional 
extrema, we first solve the mean-field equations self- 
consistently to find ip(X) solutions at f](X) = 0. We 
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seek solutions with period a, where a is the distance be- 
tween the soliton centers that are located at the mid- 
points of our primitive cells. For weak tunneling the 
solitons are wide and are accurately [40] approximated 
by solving the sine-Gordon equation that is obtained by 
minimizing the gradient approximation to the 77 = en- 
ergy functional. It follows that for dilute solitons at small 
A S as, <p(X) ~ 4 tan' 1 cxp[-( X - a/2)/£ ] modulo 2tt 
within each cell, where £ = i^J 'Airp s / 'As as is the width 
of the soliton. We find numerically that as the density 
of solitons increases (the unit cell period a decreases) 
and the energy per cross-sectional length increases, we 
eventually reach a situation at a = a c for which the 
soliton lattice solution is no longer metastable and the 
self-consistency procedure does not converge. This point 
is most conveniently identified numerically by reducing a 
for a given As^s until no mean-field solutions with finite 
phase winding exist; we discuss a more systematic but nu- 
merically more cumbersome approach below. The result 
of this calculation is summarized in Fig. 12 in which the 
regime of stable staggered-current-carrying soliton lattice 
states is plotted as a function of Ag^s and the soliton 
spacing a for d = I. For As as 0, the soliton width di- 
verges and the phase <p(X) is described by a plane-wave 
with wavevector Q. The phase diagram that we calculate 
here corresponds in this limit to the K w stability crite- 
rion and the maximum wavevector Q„, agrees with that 
found earlier (see Figs. 4, 6) to within numerical accuracy. 

0.03 1 ■ 1 ■ 1 ■ 1 ■ 1 



which metastable extrema of the energy functional ex- 
ist. The physics behind this result is not difficult to 
understand. For weak tunneling, the phase gradient, 
dx<fi ~ 27r/£ oc ^/ Asas of an isolated soliton (a ^> £) 
is below the critical value of Q* plotted as a dashed line 
in the inset of Fig. 6. In this regime the stability of 
the soliton state is therefore predominantly controlled by 
the soliton spacing a, which, for a sufficiently dense soli- 
ton lattice leads to phase gradient larger than Q*. In 
contrast, for sufficiently large interlaycr tunneling As^s 
even an isolated soliton becomes unstable, when its maxi- 
mum phase gradient dxf ~ 27r/£ oc ^Asas exceeds Q*. 
Consequently, in the large Ag^g, large soliton-spacing 
(a) regime, the phase boundary is not sensitive to soliton 
spacing a, as reflected by the approximately horizontal 
orientation illustrated in Fig. 12. We emphasize that the 
phase diagram constructed in this way does not account 
for K zz instabilities. 

The soliton lattice state carries a staggered supercur- 
rent which can be evaluated from Eq. 7 given mean-field- 
theory numerical results for the dependence of the en- 
ergy per unit length perpendicular to the current direc- 
tion e = 8/9^-/£ on soliton density. (The experimental 
significance of Eq. 7 is discussed at greater length in Sec- 
tion VI.) For each Asas, the energy per unit length for 
low soliton density (a — > 00) is proportional to the soli- 
ton density and the supercurrent therefore approaches a 
minimum value as a — ► 00: J c \ = (4/7r)^//j s Asas/k 3 - 
[48,19] Evidently J c \ decreases by decreasing the inter- 

1 11 

layer tunneling, and it vanishes like A SAS as As as — > 0. 
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FIG. 12. The stability region for soliton-lattice states, plot- 
ted as a function of Asas and soliton spacing a for d = I. This 
phase diagram is constructed by solving the mean-field equa- 
tions. Metastability in the n = plane holds in the region 
labeled Stable Soliton Lattice. 

This phase boundary shape can be summarized by say- 
ing that there is both a minimum distance between soli- 
tons and a maximum interlayer tunneling amplitude at 



FIG. 13. Maximum and minimum for the persistent cur- 
rent density of soliton lattice states in quantum Hall ferro- 
magnets vs. Asas- Ja oc ^Asas for Asas — » 0. 

At large staggered current we also find a maximum 
sustainable supercurrent density J C 2, that corresponds 
to a large soliton density. In this dense regime, defined 
by (/a > 1 solitons overlap and current is nearly uni- 
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form, leading to a highly oscillating tunneling energy that 
therefore averages away to zero. Consequently, we expect 
that J C 2 to be well-approximated by J* arising from our 
earlier criteria of K vv = for A ,5^5 = 0. Consistent 
with this, we find that J C 2, shown in Fig. 13, is not very 
sensitive to the tunneling matrix element, increasing only 
slowly with As^s. 

The experimental meaning of J c \ is unclear since one 
might expect to find solutions to generalized sine-Gordon 
equations for which the spatially averaged staggered su- 
percurrent in a finite length systems has an arbitrarily 
small value, for example by having only a partial soliton 
across the entire system. A staggered current imposed 
at the edge of the system would fix dx<px\ oc J s only 
near the boundaries of the bilayer. Provided that these 
boundary conditions are physically realistic, a question 
that would require further microscopic analysis to ver- 
ify, extrema of the functional exist for all current values 
at the edge, and for all spatially averaged staggered cur- 
rents. We discuss this point again in Section VI. 



B. K zz and K vv staggered current instabilities 

A more systematic method to determine whether or 
not a solution of the mean-field equations is metastable is 
to calculate the eigenvalues of both matrices K ZZ (X, A'), 
and K ipip (X,X'), evaluated at the corresponding energy 
functional extremum. We have made some progress with 
this approach, although we have found it to be (unsur- 
prisingly) numerically cumbersome. To make the calcu- 
lation feasible, we take the length L y of the system in 
the y direction to be finite, thereby discretizing the set 
of allowed guiding centers along the x-axis. For N guid- 
ing centers per soliton, K ZZ (X, A') and K lpip (X,X') are 
N x N matrices. In this study we increase N as far as is 
practical, then attempt to extrapolate to N = oo. Our 
numerical results are summarized in Figs. (14-17) for a 
bilayer electron system with d = I. 

In discussing these results, it is instructive to start with 
Asas = 0, since we can compare with the results ob- 
tained analytically in previous sections by Fourier trans- 
formation. In Fig. 14, we present spectrum of K zz , and 
K vv evaluated at Asas = and d — £. Starting from 
a large soliton spacing a (low staggered current), we see 
that the lowest eigenvalue of K zz decreases with decreas- 
ing a, finally vanishing for a ~ 5.2£, consistent with the 
value of Q at which the K zz instability occurs in Fig. 5. 
Because of the global phase invariance at Asas = 0, 
there should always be one zero eigenvalue of K vv at 
any value of a, an expectation that we verify within nu- 
merical accuracy. It appears from Fig. 14, that the finite 
values of L y that we use for the case of inter-soliton dis- 
tance a = 40£ (L y — 20tt£) are sufficiently short that 
we introduce artificial instabilities that do not occur for 



L y — > oo. For smaller values of a, we are able to con- 
sider larger values of L y without making our matrices 
inconveniently large. These calculations were performed 
with L y — 800ir£ 2 /a. For a = 40£ we have verified that 
the negative eigenvalue becomes smaller in magnitude, 
decreasing toward zero as L y is increased. Consistent 
with the analytical results, we find a large number of 
higher eigenvalues of K V(p approach zero as a approaches 
a* ~ 3.91, the value of a = 2n/Q at which the phase 
instability occurs for Ag^g = 0. The smallest magni- 
tude non-zero eigenvalue of K V(p becomes negative at 
a* = 3.9£, in reasonable agreement with the location of 
the supercurrent peak in Fig. 4. 
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FIG. 14. Stability matrix eigenvalues for N = 400 and 
Asas = as a function of phase period a. The pairing 
wavevector Q = 2%/a. 

In Fig. 15 we show how the stability matrix spectra 
depend on As as- We find that for a given value of a, 
increasing Asas always leads to an instability. For 
this finding is consistent with the mean-field calculation 
results discussed above. The present results indicate that 
K zz becomes unstable near the same value of Asas as 
Ktp,tp] because of the finite values we must use for L y we 
have not been able to determine which instability occurs 
first with a high degree of certainty. Figs. 15-17, indi- 
cate that negative eigenvalues in K zz do appear quite 
generally when the tunneling energy increases. The in- 
stabilities seem to appear first in K zz for small a and 
first in K VtV for large a. At a = 10£ the smallest eigen- 
value of K zz becomes negative for Asas ~ 0.028e 2 /(d?), 
near where our the mean-field theory Ly — > oo results 
indicate the first K ViV instability. Interestingly, these 
instabilities at finite As^s have no counterpart in the 
sine-Gordon theory of the soliton lattice state. They 
can be understood as being analogous to the maximum 
currents and pairing wavevectors Q that we find in the 
translationally invariant case. As Asas increases the 
maximum rate of phase winding at the center of a soli- 
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ton, Qeff = dtp/dX, increases. For an isolated soliton 
in the sine-Gordon model Q e ff ~ 27r/^ oc \/Asas- Thus 
a maximum tunneling amplitude Asas ~ 0.3e 2 /e£ is not 
unexpected given our results for the homogeneous case. 
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FIG. 15. Stability matrix eigenvalues for N = 400 and 
a — 101 as a function of Asas- 
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FIG. 17. Stability matrix eigenvalues for N = 400, and 
a — 40£ as a function of Asas- 



VI. DISCUSSION AND CONCLUSIONS 
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FIG. 16. Stability matrix eigenvalues for N = 400 and 
a — 20£ as a function of Asas- 



In this paper we have presented a mean-field theory of 
staggered supercurrent carrying states in quantum Hall 
bilayer superfluids. Our analysis was based primarily on 
the behavior of two stiffnesses K zz and K vv that char- 
acterize respectively the energy cost of small interlayer 
charge and interlayer phase fluctuations. By comput- 
ing these quantities, we have been able to address the 
physics that controls the dependence of critical staggered 
current on interlayer tunneling amplitude, layer spacing, 
and other tunable parameters that characterize bilayer 
systems. We also studied the characteristics of weak links 
in quantum Hall superfluids that are created by a gate 
voltage. Unlike the field-theoretic approach [19], the mi- 
croscopic analysis here does not rely on a small wavevec- 
tor expansion, and therefore has the virtue of accounting 
for the full wavevector dependence of the stability matri- 
ces. This is especially important because the wavevector 
expansion of these matrices is not analytic which indi- 
rectly causes these quantities to have minima at finite 
wavevectors. The staggered critical supercurrent predic- 
tions presented in this paper are consequently more real- 
istic than those derived from the field-theoretic approach 
[19]. 
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FIG. 18. Schematic representation of the staggered cur- 
rent-induced QH supersolid state for which periodic interlayer 
charge modulation and QH superfluidity coexist. The peri- 
odic rate of phase winding is also illustrated. 

There are a number of limitations to the approach we 
have followed. Our analysis is confined to harmonic in- 
stabilities and will therefore miss any that are induced 
through a nonlinear mode, nucleation of vortices for ex- 
ample. In addition, our harmonic analysis does not un- 
ambiguously determine the state to which the system 
evolves when the instability point is reached. Based 
on experience with ordinary superfluids and supercon- 
ductors, a natural guess is that once one of the stiff- 
nesses vanishes, the current-carrying superfluid state of 
the quantum Hall bilayer is unstable to a resistive (nor- 
mal) state in which staggered current leads to a stag- 
gered voltage. We believe that this scenario applies when 
K w controls the instability, and when the instability in 
K zz (q) takes place at a large wavevector corresponding 
to a length scale of order or smaller than the magnetic 
length I. The instability process has that character for 
the large interlayer spacing values, d > £, relevant to cur- 
rent experiments. However, in the opposite limit, d << £, 
it is likely (based on similarity with our earlier work on 
the interlayer-charge imbalance instability driven by an 
in-plane magnetic field [50,37]) that vanishing K zz (q*) 
instead signals an interlayer charge imbalance quantum 
(nonequilibrium) transition to a distinct staggered cur- 
rent carrying QH superfluid state. As with the equilib- 
rium canting transition [50,37] at small d/£, the instabil- 
ity signaled by vanishing of K zz (q*) occurs at a small but 
finite wavevector q* (with q*i « 1), and therefore cor- 



The vector potential A has no physical meaning if 
Asas = 0. For A S as 7^ 0, however, an in-plane field 
is represented by A = B\\d and in equilibrium Eq. 14 de- 



responds to a development of a dipole-stripe state akin to 
a unidirectional charge density state with period 27r/g*, 
as illustrated in Fig. 18. Since, in addition to staggered 
gauge symmetry, such a state also spontaneous breaks 
translational invariance along the staggered current di- 
rection, it corresponds to a quantum Hall supersolid with 
coexisting superfluid, solid, and quantum Hall orders. 



J,(X), a 
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FIG. 19. Schematic QH bilayer in the presence of finite sin- 
gle-particle interlayer tunneling Ag^s- When there is collec- 
tive tunneling between the layers, the staggered supercurrent 
must be spatially dependent in order to conserve the difference 
in charge between the layers. The spatially averaged stag- 
gered supercurrent can be arbitrarily small in a finite systems 
becaus e of edge e ffects. The soliton width in this illustration 
8 = iyj p s / Asas, while the period a is the distance between 
soliton centers. This schematic plot shows a partial soliton 
on the right that can make an arbitrarily small contribution 
to the spatially averaged staggered supercurrent. 

Another feature of our work is the finding that 
the staggered super current evaluated for soliton-lattice 
states using Eq. 7 has a minimum value J c i as well as a 
maximum value. This expression for the staggered cur- 
rent can be motivated in several different ways, but its 
connection to experiment is not always direct, as we now 
explain. The general expression for the staggered current 
can be derived by introducing a spatially constant stag- 
gered vector potential = — A± = A/2 and noting that 
the interaction energy is then a function of the gauge 
invariant inter-layer phase difference fx + eXA/Tic. It 
follows that the staggered current is given by 



(14) 



scribes the persistent currents responsible for orbital dia- 
magnetism. In the absence of an in-plane field A = 0, and 
<p{X) is determined by minimizing Eq. 2. When the mini- 
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14 



mization procedure is carried out with periodic boundary 
conditions, as in our numerical calculations, changing A 
is equivalent to changing the period over which <p is re- 
quired to change by 2ir, and Eq. 7 follows. For realistic 
systems, however, edge effects can have a large effect on 
the current, especially when individual solitons, which 
carry circulating currents in their cores, are well sepa- 
rated. 

J c i is analogous to the lower critical magnetic field 
H c i in typc-11 superconductors. Only for J a > J c \ does 
it pay for the system to give up some of the tunneling 
kinetic energy and allow a soliton that carries this mini- 
mum staggered current J c \ into the system. J c \ is closely 
related to the critical in-plane magnetic field that induces 
the commensurate-incommensurate transitions in quan- 
tum Hall bilaycrs. [9,50] 

Consider, for example, the limit of slowly varying 
phases for which the local staggered current density 
can be shown to be [9] 2p(dip/dx), in agreement with 
Eq. 14. In this case the spatially averaged staggered 
current is proportioanl to the number of solitons in the 
system, which can clearly be arbitrarily small. This 
conclusion should be contrasted with that reached by 
Shevchenko [48], working with a closely related long- 
wavelength model that applies to the case of of electron- 
hole bilayers at zero magnetic field. (Excitonic super- 
fluidity has not yet been conclusively established in ex- 
periment for this case.) Shevchcnko's conclusion that 
J c i would represent an experimental minimum staggered 
supercurrent was based on the argument that finite cur- 
rent configurations could be studied by adding a soliton 
chemical potential term [more precisely a "chemical po- 
tential" for interlayer phase twist] J s 2p s dx<f/fo to the 
bilayer energy density functional. This term is chosen so 
that the modified energy functional has its absolute min- 
imum when Eq. 7 for the current density is satisfied, and 
is modelled after related approximations that are com- 
mon [49] in the description of current-biased Josephson 
junctions. In this approach, because for finite tunnel- 
ing As^s the induced phase twist costs interlayer kinetic 
energy (tunneling), a minimum staggered current J c \ is 
required for phase winding to occur. Below this critical 
value, J s < J c i (analogous to a Meissner state of type- 
II superconductor, where vortex density cx B vanishes) 
no solitons are induced, and phase ipx remains spatially 
uniform and pinned at (except for the boundaries). 
The consequent vanishing of the interlayer tunneling cur- 
rent leads to a conserved staggered current J s (x), that 
is therefore spatially uniform. In our view, however, this 
thermodynamic description of finite current states does 
not apply to electron-hole bilayer systems with excitonic 
superfluidity, cither at zero field or in the quantum Hall 
regime. The key difference between quantum Hall bilay- 
ers and Josephson junctions is the absence of a supercon- 
ducting current-bias environment which imposes a defi- 
nite value of the staggered curent at every point in space. 



In the bilayer quantum Hall case, current-carrying states 
are never equilibrium states, but they can be metastable. 
External contacts can in principle, impose a given value 
of the staggered current at every point on the edge of 
the system, assuming that non-collective current carry- 
ing excitations can be completely discounted. However, 
the staggered current density at interior points must be 
determined by finding local minima of the energy func- 
tional that are consistent with edge boundary conditions. 
In general when As^s is non-zero, many such solutions 
exist. The critical current is naturally defined as the 
metastable configuration with the largest staggered su- 
percurrent, even for geometries that are more general 
than those considered here. 

It is difficult for several reasons to use our microscopic 
HFA approach to calculate the activation energies for the 
saddle point extrema that separate local minima of the 
HF energy with different numbers of phase slips. These 
activation energies will control the thermal phase-slip nu- 
cleation rate and the associated staggered voltage at any 
finite temperature. One problem encountered in seek- 
ing these solutions is that standard simple iterative ap- 
proaches can be used to find solutions of the mean-field 
equations that correspond to local minima, whereas some 
more subtle and clever technique would have to be used 
to coax the system iteratively toward a saddle-point so- 
lution. A second problem is that for the wide Hall bars of 
interest to us here, the saddle point solution likely con- 
sists of nucleating a vortex and moving it to the middle of 
the Hall bar; these solutions have polar and azimuthal an- 
gle fields that depend on both spatial coordinates, mak- 
ing the problem two-dimensional. 

We can, however, make some progress toward under- 
standing staggered channel dissipation by combining our 
HFA results and a scaling analysis supported by exact 
field-theoretic calculations valid for a short-range inter- 
action model [19]. For example, consider the case of 
thermally-activated staggered resistance controlled by a 
saddle-point phase-slip solution m s z p (x), similar to the 
one illustrated in Fig. 20. We expect that the limiting 
saddle point will have this character for narrow Hall bars. 
The phase slip is characterized by a single length scale 

v 1 - ml 



which reduces to the width £ = \J Ps/K zz (0) [19] of 
saddle-point solution in the limit of short-range interac- 
tions or small £09*7 valid for d/l <C 1; q* is the wave vector 
at which, in the presence of long-range Coulomb interac- 
tion K zz (q) is minimum. In this short-range limit, we ex- 
pect the saddle-point solution mf (x) and the associated 
energy barrier to be well described by the field-theoretic 
prediction [19]. 
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FIG. 20. Schematic illustration of the saddle-point solu- 
tion, m s z p (x) for £og„ — » 1~, characterized by envelope width 
£ and wavevector . The free-energy associated with this 
m s /(x) controls the energy barriers for thermal activation of 
phase slips and dissipation in the staggered current channel 
when the Hall bar is narrow. 

In the opposite limit, as the width £ of the soliton nu- 
cleation site approaches (from below) the period at 
which the intcrlayer charge-imbalance m z = tit — ub in- 
stability takes place, we expect a saddle-point solution to 
be qualitatively modified by long-range Coulomb inter- 
actions. In this regime we predict an oscillatory m B z p (x), 
with wavevector q* and an envelope of width £, as illus- 
trated in Fig. 20. The associated energy barrier in the 
limit of vanishing staggered current is given by 

E B =el{l-$ql) 1/2 L yi (16) 

where e% the energy barrier per unit of soliton length in 
the limit of short-range interactions [19]. 

This analysis has been limited so far to one- 
dimensional instabilities. The energy barrier associated 
with such a one-dimensional soliton grows with the width 
L y of the quantum Hall bar. In contrast, for vanishing 
tunneling A$as = 0, barriers for pomi-vortex nucleation 
are proportional to p Sl up to corrections logarithmic in 
L y . Consequently, for sufficiently wide quantum Hall 
bars, L y > £ [19], the staggered resistivity will be deter- 
mined by point-vortex nucleation and unbinding. Indeed 
standard arguments guarantee that some dissipation will 
always be present at finite temperatures and voltages, as 
we discuss below. 

Our discussion has so far ignored quantum and ther- 
mal fluctuations and quenched disorder, present in real 
systems. For layer separations ~ 10% below the criti- 
cal value d c [52,53], quantum fluctuations in the absence 
of quenched disorder can be treated perturbatively [52] 
and will yield only small corrections to our results. In 
contrast, thermal fluctuations and especially in combina- 



tion with disorder can modify some of our results qual- 
itatively. Even in the absence of quenched disorder, for 
vanishing tunneling, thermal fluctuations lead to a finite, 
activated dissipation (i.e., staggered voltage) at any fi- 
nite staggered current, and therefore as in the case of 
superconductors, preclude an unambiguous definition of 
a critical current. In this case a supcrfluid state is distin- 
guished from a dissipative one by a vanishing linear resis- 
tivity. Extending standard arguments [51] to our system 
for Ag^s = limit, we predict a nonlinear power-law 
staggered I-V characteristic E s <~ Jf, a(T) > 3 to char- 
acterize the quantum Hall staggered superfluid transport. 
Despite finite dissipation even below mean-field critical 
staggered current J c computed here, at low temperatures 
we nevertheless expect a strong crossover in the I-V at J c , 
with the staggered voltage drastically dropping for cur- 
rent below J c and dissipation controlled by thermally ac- 
tivated vortex nucleation and vortex transport over bar- 
riers. 




FIG. 21. Low temperature and weak disorder staggered 
current-voltage (I-V) characteristics for QH bilayers for (a) 
finite and (b) vanishing interlayer tunneling Asas- 
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For finite interlayer tunneling, As^s 7^ 0, the ion- 
ization energy necessary to separate vortex pairs grows 
linearly with separation R as (psAsas) 1 ^ 2 R/ Con- 
sequently, in thermodynamic limit (wide bilayers) wc 
should expect a true staggered critical current J c \ = 
(l/M)(p s AsAs) 1 ^ 2 7 below which staggered transport is 
dissipationless, E(J S < J c i) = 0, even at finite tempera- 
ture, and for J s > J cl , E(J) oc \J S - J c i| Q ( T ' [19]. The 
corresponding I-V is illustrated in Fig. 21a. 

Quenched disorder is another important ingredient 
that we have not considered in detail here. As we have al- 
ready emphasized, our analysis of the weak-link strongly 
suggests that by creating weak links, disorder will sup- 
press critical currents. Short-range disorder potential 
fluctuations will tend to nucleate vortices, a fact that 
can be understood in terms of the short-range-correlated 
random effective vector potential [54,25] that it gives rise 
to. For Asas — * 0, the random vector potential model is 
known to have a stable low temperature quasi-long-range 
ordered phase with no unbound vortices that appears be- 
low a critical strength of disorder. [55] (Microscopic cal- 
culations [52] that treat disorder using a self-consistent 
Born approximation, lead to the same conclusion.) Be- 
yond the presumed disorder threshold, the system is in 
the gauge-glass regime, which in two dimensions is un- 
stable to vortices that are mobile at finite temperature 
and only localize as T — > 0. Hence in the weak disorder 
limit, when vortices are dilute, <5 <C Ll (where Ll is the 
Larkin length), the mean- field critical current found in 
this paper should be clearly observable. In fact for weak 
disorder we expect to see three regimes of dissipation 
via motion of disorder-induced vortices: (i) slow, ther- 
mally activated creep for J s < J p , (ii) vortex flow, for 
J p < J s < J c (where J c is the mean-field critical current 
of a clean system derived here) , and (iii) normal state 
dissipation for J c < J s , when the system is past its point 
of metastability. The corresponding I-V is schematically 
illustrated in Fig. 21b. 

How then might the critical currents calculated here, 
which are based on a Landau-like metastability require- 
ment for the staggered-current carrying state, appear 
experimentally? The best possibility appears to be by 
looking for features in the non-linear I-V relationships 
seen in interlayer quantum Hall drag experiments. [23] 
So far these experiments seem to show a linear relation- 
ship between the longitudinal staggered voltage and the 
staggered current, with a thermally activated dissipation 
process. There is no hint at present that the staggered 
electric field initially grows sublincarly with current, sug- 
gesting that J p might be driven to zero by disorder. It 
will be very interesting to see whether or not a relatively 
sharp crossover can be observed experimentally at low 
temperatures at a current density that can be associated 
with J c . The crossover would be from a regime in which 
the staggered current is largely collective, but not com- 
pletely dissipationless due to the motion of vortices nu- 



cleated thermally and by disorder, to a regime in which 
there is essentially no collective staggered current. If our 
estimates are qualitatively reliable, the scale of J c is nor- 
mally larger than those used in typical experiments but 
vanishes as (d c — d) 1 ' 2 as the phase boundary for the 
spontaneous coherence state is approached from the be- 
low. For this reason, we believe that nonlinear transport 
experiments in the drag geometry, current in one layer 
and voltages measured in both layers, as a function of 
d c — d will be helpful in sorting out the complex physics 
present in current samples. [23] 
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